Rapid purification of quantum systems by measuring in a feedback-controlled 

unbiased basis 
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Rapid-purification by feedback — specifically, reducing the mean impurity faster than by mea- 
surement alone — can be achieved by choosing the eigenbasis of the density matrix to be unbiased 
relative to the measurement basis. Here we further examine the protocol introduced by Combes and 
Jacobs [Phys. Rev. Lett. 96, 010504 (2006)] involving continuous measurement of the observable J^ 
for a D-dimensional system. We rigorously re-derive the lower bound (2/3) (D -f 1) on the achiev- 
able speed-up factor, and also an upper bound, namely -D^/2, for all feedback protocols that use 
measurements in unbiassed bases. Finally we extend our results to n independent measurements 
on a register of n qubits, and derive an upper bound on the achievable speed-up factor that scales 
linearly with n. 



PACS numbers: 03.67.-a,02.30.Yy,02.50.-r,89.70.-a 



I. INTRODUCTION 

Many quantum information processing (QIP) tasks re- 
quire pure states as a resource [l|. Preparation of pure 
states may be a bottleneck in certain physical implemen- 
tations of some QIP protocols {e.g. solid state systems). 
In these systems measurement can not be treated as in- 
stantaneous 0Tal' ^^^ ^1'*^ often slow compared to the 
timescale for unitary operations [6[. In such systems, a 
projective measurement is approached only in the limit 
of long measurement times, by integrating the measured 
current Q- There are demonstrable advantages associ- 
ated with using the dynamics of the measured current 
to perform read-out or state estimation Q. Technology 
is approaching the point where it can implement useful 
feedback protocols [8l-[la|. Thus it makes sense to con- 
sider combining continuous measurement with quantum 
feedback [1^ to control the purification process. Our goal 
is to increase the rate at which these systems can be puri- 
fied in order to alleviate the bottleneck identified above. 
We note that, by the addition of a unitary at the end of 
the measurement process, ra pid purification is identical 
with rapid state preparation |17H19j . 

The Von-Neumann entropy, S{p) = — Tr[plogp], for 
a quantum state p is a natural choice for characteriz- 
ing a procedure that purifies and/or cools the state of 
a quantum system. Unfortunately, the von Neumann 
entropy is not an easy quantity with which to perform 
calculations. Because of this the linear entropy, given by 
L[t) = 1 — Tr [p(t)^] , is often used to make the calcu- 
lations more tractable. In this paper we use the linear 
entropy, also known as the impurity, to characterize the 
purification process. Since unitary operations on a quan- 
tum system leave the eigenvalues of the density matrix 
unchanged, in order to reduce the entropy of a system it 
must be allowed to interact with a bath. It is this inter- 
action which enables, and is common to, all purification 



and cooling procedures. 

Until recently the majority of work in the field of state- 
purification had been based on what we will term 'open- 
loop' control methods, such as algorithmic cooling [20| or 
dissipation engineering [2l|. In this approach the bath 
remains unmeasured or the results of the measurements 
arc forgotten. As a result the system evolves according 
to a deterministic master equation, such as 

dp = dtC [P\ p = -idt[H{t),p] -t- 27 dtV[e] p, (1) 

where V[A]p= ApA^ - \{A^ Ap + pA^ A) [l|, H{t) is 
the system Hamiltonian, I is a, lowering operator and 7 
is a relaxation rate. Consequently the control signals or 
interactions are not conditioned on measurement results, 
and thus can be completely determined prior to the cool- 
ing process. Cooling can be thought of as a version of 
purification in which the final state is also the ground 
state for the system. 

In measurement-based purification schemes the bath 
must be measured. In this case the system can be puri- 
fied without control fields, provided the system observ- 
able coupled to the bath has no degenerate eigenvalues, 
by measurement alone. The change to an observer's state 
of knowledge of an individual system after a weak mea- 
surement of the system observable X is described by the 
stochastic master equation (SME) [22, [23| 



dp = dtCc [X] p = 2jdtV[X]p+ v/27 dW H [X] p, (2) 

where H[A]p = Ap + pA^ - Tr [{A^ + A)p] p [l6l. and 
dW is the increment of a Wiener noise process [2J]. It 
should be noted that we have moved to a frame that has 
enabled us to factor out the Hamiltonian evolution. The 
measurement strength, 7, determines the rate at which 
information is extracted, and thus the rate at which the 
system is projected onto a single eigenstate of X [23, US, ■ 
(This means for times r ^ 7"^ we may say that wc have 



performed a projective measurement of the observable 
X.) The measurement resuh in a smah time interval 

[t, t + dt) is 



(3) 



dR = ^A'y{X{t))dt + dW{t), 



where dW is the same Wiener noise process that appears 
in Eq.©, and {X{t)) = Tr[Xp(i)]. We wih denote the 
continuous measurement record obtained by the observer 
integrated up until time t as R{t) = J^ dR{t'). 

Conceptually a continuous measurement is a sequence 
of weak measurements, in the limit where the strength of 
the measurements tends to zero and the repetition rate 
tends to infinity [H [H H, [l^j. Feedback, if thought 
about in the framework of a sequence of weak measure- 
ments, amounts to the ability to perform a unitary oper- 
ation in-between each measurement [23l- While unitary 
operations by themselves do not change the entropy of a 
system, the average reduction in the linear entropy caused 
by a measurement depends not only on the measurement 
but also on the prior state of the system. Feedback proto- 
cols that increase the rate of purification usually involve 
choosing the applied unitary after each weak measure- 
ment so as to increase the entropy reduction generated 
by the next weak measurement [3, I17H191 128| - |43| . 

Here we consider rapid purification feedback protocols 
that apply unitary operations to continually keep the 
eigenbasis of the density matrix unbiased with respect 
to the basis of the observable being measured, X. (Two 
bases of a D-dimensional system are unbiased with re- 
spect to each other if the inner product of every vector 
of the first basis with every vector of the second basis 
is equal to l/VD.) We will refer to this class of feed- 
back protocols as "unbiased-basis" (UBB) feedback pro- 
tocols. Specifically, the goal of these protocols is to in- 
crease the rate at which the average linear entropy of 
the system, {L{t)), decreases as a function of time. The 
average here is taken over all possible trajectories (mea- 
surement records) for the evolution of the system under 
the measurement (and feedback). For a single qubit, the 
UBB protocol is optimal [H [H El- UBB protocols 
were first examined for systems of arbitrary dimension 
in [31| . The purpose of the present work is to clarify and 
extend the results in [3l|, as well as to derive upper and 
lower bounds on the performance of UBB protocols when 
applied to registers of qubits. 

In Section |lT] we review rapid purification for a sin- 
gle qubit. Section [llll discusses previous results on qudit 
feedback in an unbiased-basis. And Section ITVl provides 
a unified formalism to derive both the upper and lower 
bound on purification for any UBB feedback on a qudit. 
Finally, we derive upper and lower bounds for UBB pro- 
tocols applied to a register of qubits in Section |Vl and 
conclude with a discussion of the results. 



II. RAPID PURIFICATION FOR A SINGLE 
QUBIT 

In this section we review the UBB protocol for a single 
qubit (originally presented in [30|). For a single qubit 
the UBB protocol is optimal - a non-rigorous proof of 
this is given in |30[, and a rigorous proof is given in 18 1 
and confirmed in Ref. |43j . It is not presently known 



whether UBB protocols are optimal for systems of higher 
dimension. 

To calculate the factor by which a feedback protocol 
speeds up the purification of a system we divide the time 
it takes the feedback protocol to achieve a given value of 
{L(t)) by the time it takes a measurement in the absence 
of feedback to achieve this value [43 • In both cases we 
start the system in the maximally mixed state. We will 
refer to a measurement acting without feedback as a bare 
measurement 44 1 . 



A. Purification from Measurement Alone 

We will perform our analysis keeping the dimension of 
the system arbitrary, for later convenience. Consider a 
quantum system of dimension D initially in the maxi- 
mally mixed state p(0) = l/D, where I is the D x D 
identity matrix. It is possible to obtain a closed form 
expression for the linear entropy as a function of the 
measurement record by using the linear trajectory for- 
mulation of continuous measurements [la, |4^, |4y| (for a 
simple introduction see [22|). This involves solving the 
linear version of the SME, which produces an unnormal- 
ized density matrix. The linear version of the SME in 
Eq. © is 



dp = dtCc [X] p = 2-fdtV[X]p+ V27 dR H [X] p. (4) 

Here H[A] p = Ap + pA^ , and the tilde over p denotes the 
lack of normalization at all but the initial time. Because 
the initial state is p = p{0) = l/D, the density matrix 
commutes with X at all times, and this makes obtaining 
the solution simple. This solution is [22 1 



p{R,t) = ciip{-A-fXH)cxp{2y^2jXR{t))I/D, (5) 

where, as before, R{t) is the integrated measurement 
record. For a qubit, for which D ~ 2, wc take the ob- 
servable to be X = J^ = 0-^/2. The solution becomes 



p{R,t) 



= -7* 









V27fl 



(6) 



(Here we have dropped the time-dependence of R for 
compactness.) 

The final normalized density matrix p{R, t) is given by 
dividing p{R, t) by its norm, M — Tr[/5(i?, t)]. The proba- 
bility density that we obtain the state p{R, t) at time t, is 
given by V{R,t) = TVexp (-i?2/2f)/757rt. The average 



impurity of the final state is thus given by averaging the 
impurity L[p{R,t)] over the probability density V{R,t). 
For a single qubit this gives 



im) 



.-7* 



^+00 



-R^/2t 



cosh(V2^i?) 



dR. 



(7) 



While this integral has no analytic solution (to our knowl- 
edge), we can obtain the behavior in the long-time limit 
by noting that the integral contains two multiplied distri- 
butions. The distribution in the numerator is broad com- 
pared to the distribution in the numerator for t 3> J~^- 
Thus in this long time (LT) limit, the integral can be ap- 
proximated by J_'^ dR/ cosh{^/2^R) = 7r/y27, and we 
have 



(Lit)) 



-ft 



LT 



y/T&TTjt 



(8) 



The key result is that the impurity for a bare continuous 
measurement scales asymptotically as e~'^*. 

We note that recently Jordan and Korotkov [3| gener- 
alized Eq. (O for an arbitrary initial state, po = p(0) = 
i(I -I- xax + ycTy + za^). Using linear trajectory theory 
as above one can show that the impurity decays as 



im) 



--^'Lipo] 



+CX3 



-«'/2*di? 



, cosh(\/27i?) + zo sinh(V27-R) ' 

. ^^) 
where zq = Tr [cr2p(0)]. The asymptotic expression is 



(i(i)>LT = 



Tre-^* 2L(po) 



vTB^yr^^ 



(10) 



B. Purification Using Feedback 

For the initial state p(0) — 1/2, the measurement dy- 
namics is symmetric with respect to rotations about the 
z-axis. Because of this we lose nothing by restricting our 
feedback protocol to rotations about the y axis at some 
rate a{t). The SME for this situation is 



dp 



-idt[a{t)Jy, p] + dtCc [Jz] pit). 



(11) 



To simplify the calculations we assume that a{t) can be 
arbitrarily large compared to the measurement strength 
7. This means that we can consider the action of the 
control in each infinitesimal time- interval, [t,t + dt), as 
a unitary that generates a rotation through any desired 
angle. This unitary is therefore of the form 



Ut = U{t + dt, t) = exp {-ia{t)Jydt}. 



(12) 



In this case the state following an infinitesimal time-step 
consisting of measurement and feedback is 



Pfh{t + dt) 



Ut{p{t) 
Ut{pit) 



dp{t)}Ul 
dtCc[J.]p{t)}Ul 



(13) 



Note that up to a unitary transformation, which has no 
effect on the purity, this is equivalent to having the feed- 
back change the measurement basis: 

Pfb(t + dt) = p{t) + dtCc[X{t)\p{t), (14) 

where X{t + dt) = UtX{t)U^ and X{0) = J^. This is 
a Heisenberg picture with respect to the control unitary. 
In what follows we will always assume, for the sake of 
simplicity, that the feedback changes the measurement 
basis, rather than the state of the system. 

To derive an expression for {L(t)), we begin by exam- 
ining the first-order change in the linear entropy. This 
is 



dL = d[l- Tr [p2]] = -Tr [d{p^)] 
= -Tr [2pdp + (dp)2] . 



(15) 



Here we must keep the second-order term in dp, because 
p{t) is stochastic, and (dW)'^ = dt ^. From Eq. ©, 



with X replaced by X, the change in impurity is thus 

dL = -87<it{Tr [pXpX] - 2Tr [Xp] Tr [Xp^] 
+ Tr[pXfTT[p'-]} 
-4^/2JdW{^T [Xp^] - Tr [pX] Tr [p^]}.(16) 

To obtain the greatest decrease in impurity in each time 
interval [t,t + dt), we must now optimize over all uni- 
taries Ut to obtain the locally optimal X{t). It was shown 
in [30| that for a single qubit this is achieved by choosing 
the eigenbasis of X to be unbiased with respect to the 
eigenbasis of p. Each infinitesimal measurement disturbs 
this unbiased relationship, and thus feedback is required 
to maintain it. Since the disturbance to the basis of p is 
proportional to dW (rather than dt), to keep the bases 
perfectly unbiased requires that a{t) be arbitrarily large. 

If the bases of X{t) and p{t) are unbiased, then the 
bases of X(t) and p{t + dt) will only be infinitesimally bi- 
ased, so that the necessary feedback unitary U{t,t + dt) 
will be infinitesimally different from I. This is required 
for physically reasonable feedback. Mathematically, how- 
ever, it is simpler to imagine the case where p{t) is diag- 
onal in the J^ basis. In this case, X{t) will be obtainable 
from Jz by a finite unitary rotation. 

We choose X = Jx,so that the unitary that transforms 
X toX viaX = TJzT'<, is 



T = exp 



(■f 



J, 



V2I-I 



(17) 



One can show that if X is traceless, and unbiased with 
respect to p, then regardless of the dimension Tr [Xp] = 
and Tr [^p^] = (the derivation is given in Ap- 
pendix [^ . This considerably simplifies the expression 
for dL, which is now 



dL = -87Tr [XpXp] dt. 



(18) 



For a single qubit it turns out that Tr [XpXp] = L/4 
and we obtain the very simple equation 



dL = -2'-fdtL. 



(19) 



Evolution of the impurity is thus deterministic, and is 
given by 



L{t) = e^27ti(o). 



(20) 



Note that the evolution of L is only perfectly determin- 
istic under the assumption that the observable X is per- 
fectly unbiased with respect to the density matrix. In the 
long-time limit we obtain the speedup factor from Eqs. 
(|5]) and (PO)) . Denoting the time taken by a bare mea- 
surement to reach a given value of L as tbaro, and that for 
the feedback protocol as ifb, we equate L(ifb) — i(ibarc) 
and solve for the ratio ifb/ibaro- Doing so gives 



^barc 



1 



In \/167r7tba 

27ibare 



ln27r 

27ibarc 



(21) 



For sufficiently large ibaro (equivalently, a sufficiently 
small target impurity) the second and third terms are 
insignificant, and we obtain 



tfh 



(22) 



Numerical calculations show that for shorter times 
(higher target impurities) the speedup factor is always 
less than this. Thus the largest possible speedup for a 
single qubit is a factor of 2. 



III. QUDIT PURIFICATION USING 
FEEDBACK 

A. Purification from Measurement Alone 

To analyze UBB feedback protocols for Z3-dimensional 
systems, we will need the asymptotic evolution of the 
impurity for a bare measurement of Jz for arbitrary D. 
In [3l| it was shown only that at long times L ^ e~^*'. 
A more detailed analysis, to be presented elsewhere [47| . 
shows that the asymptotic behavior is 



{L2{t))uT = 



2{D-l) 



.-It 



D 



^IQ^tlT 



(23) 



Interestingly this is merely the expression for a single 
qubit, Eq. ([7]), multiplied by a factor that depends on 
the dimension of the system. 



B. The UBB protocol 

In this section we examine the generalization of the 
qubit rapid purification algorithm to D-dimesional sys- 
tems (qudits), which was first proposed in |3l| . Recall 



that an UBB protocol for a D dimensional system in- 
volves using feedback to continually adjust the measured 
observable X(t) so that its eigenbasis remains unbiased 
with respect to the density matrix. Note that for these 
protocols we require only a single basis that is unbiased 
with respect to p — we do not require a complete set of 
mutually unbiased bases. 

The change in the impurity in a single infinitesimal 
time-step is once again given by Eq. (|18p : 

(24) 



dL = -87Tr [XpXp] dt. 



However, in this case the right-hand side is no longer a 
simple function of L. Further, the evolution of L need 
no longer be deterministic; L is now coupled to other 
functions of p, and these will still evolve stochastically. 
It will be useful in what follows to write the right-hand 
side explicitly in terms of the matrix elements of X in 
the elementary basis (the eigenbasis of p) . Denoting this 
basis as {\i)}, i = 1,. . . ,N, we have 



dL = -87dt^|X,,jf A,Aj. 



(25) 



where Xj is the eigenvalue of p associated with the eigcn- 
state \j). 

In 31| a lower bound was derived on the performance 
of UBB protocols, by considering a protocol in which 
simultaneous measurements are made of all observables 
that can be obtained from X by permuting the basis vec- 
tors. This was possible because the evolution of L under 
this protocol can be solved. Under the Dl simultaneous 
measurements, the total change in L, dLtot, is simply the 
sum of the dL's due to each measurement: 



dLt. 



8jdt 



D\ 



J2 Tr [xIpXLp] 



(26) 



?n — 1 



where the X,n— Pm-^-P-rn are D\ permutations of the op- 
erator X. Since at least one of the observables Xm in the 
sum in Eg. ([25)) must give a dL that is at least as large 
as the average over all the Xm, the performance of this 
protocol is a lower bound on the performance of proto- 
cols that employ a single optimized observable X{t). An 
intricate calculation j31[ shows that Eq. ([25]) can then be 
rewritten in a remarkably simple way: 

dLtot = -S-fdt^^^^p^Ti [X^] L. (27) 

In Appendix [Bj we give a new, and much more detailed, 
proof of the process to obtain the relation in Eq. ([27)) from 
Eq. ([26]). and in Sec.|V]we will use this to obtain a lower 
bound on UBB protocols for a register of qubits. When 
X = Jz, Tr [^2] = Tr [j^] = D{D^~1)/12, and Eq. ^ 
becomes 



dL 



LB 



ht{D + l)L. 



(28) 



The subscript LB indicates that this increment is a lower 
bound on \dL\ for any unbiased basis feedback, that is 
\dLi^B\ < \dL[p]\. 



IV. UPPER AND LOWER BOUNDS FOR 

RAPID PURIFICATION USING UBB 

FEEDBACK 

We now introduce a new method to determine the 
lower bomid obtained in [31|, Eq. (|28p here, and a similar 
method that allows us to obtain an upper bound. 

The lower bound in Ref.|3l| was obtained by averag- 
ing over all permutations Pm of the measured observable. 
We now note that this procedure renders (dL) invariant 
to such permutations. In light of this, and by analogy 
with the technique employed in Ref. |4l| , we introduce a 
density matrix for which dL is invariant under the per- 
mutations Prr, : 



(29) 



PF = diag I 1 - A, -^, -^, . . . , -^ 



We call this the "flat state" , as it has one large eigenvalue 
and the remaining eigenvalues are equal in magnitude 
( "flat" ). For any Z?-dimensional state p, with an impurity 
L[p], we can always find a A so that ^[pf] = L[p\. 

Intuitively the state in which dL is most sensitive to 
permutations for a fixed value of L is 



P2 = diag(l-A',A',0,...,0). 



(30) 



We will refer to this state as the "binary" distribution. 
Once again, for a given p we can always find a A' so that 
L[p2]=L[p]. ^ 

We now derive upper and lower bounds on UBB pro- 
tocols by showing that when L[p] ~ L[p-p] ~ L[p2], and 
when each of these three density matrices have a permu- 
tation applied that maximizes dL, 



\{dL[pT.])\<\{dL[p])\<\{dL[p2 



(31) 



The Lower Bound via the flat distribution 



In Sec. mil we showed that the lower bound on the 
change in impurity is dLi^B = ~^lt{D + l)L. That is, 
M-^lbI < \d'L{p)\. Now we show that dLLB = dL{pp). 
We first note that the impurity of pp can be written as 



L[pf] = {D-l) 



2A(1-A) 



D- 1 



{D 



{D-l) 



(32) 



Equation (j32[) will enable us to factor out the impurity 
in the following working. Substituting pp into Eq. (|25|) 
to calculate dL{pF) (or dLp in shorthand notation) and 
then writing dLp in the form of Eq. p2p we have 



dLt 



'8'y dt 



/ 2A(1-A) ^ 



where the largest eigenvalue 1 — A is associated with some 
particular eigenstate \p) . It is easy to see that dLp is in- 
variant under the transformation dL[pi?] — dL{P^j^ppPr,i) 
for all the D\ permutations Pm labeled by m. The per- 
mutation invariance oi pp already implies that dL{pp) = 

(l/-D!)E™Li^^(^rtPf-P™) = dLtot - d^LB. Neverthe- 
less we persue the simplification of Eq. p3p by noting 
that for any unbiased-basis 

Y^lXrpf = {D + 1){D-1)/12, (34) 

as shown in Appendix [U] From this it follows that 



7 , \^r,c 



\^^£±l{D-l){D-2). (35) 



Thus it is possible to write Eq. (|33l) so that it has the 
same form as Eq. (I32p . simplifying the resulting expres- 
sion gives: 



dLf 



~8jdt 



JD + l) 
12 



L[pp]^-liD- 



l)-fdtL{t). (36) 



It should be noted that technically we have not inde- 
pendently rederived the lower bound in Eq. (|28p: rather 
we have shown that for all unbiased bases and for all pos- 
sible permutations Pm of pf's basis and for all impurities 
dLi^B = dLp. 



B. The Upper Bound via the binary distribution 

We now show that |di[p]| < |(iL[p2]|- Substituting p2 
in to Eq. ^ gives dL2 = -8jdt[2{l - A')A'|X,cP]- 
This expression is sensitive to the arrangement of the 
eigenvalues of p2; accordingly dL[p2] ^ dL{P^^p2Pm) for 
most permutations. The impurity for p2 can be written 
as 



L[p2] = 2(1 - A)'A' = L[p] = Y, ^rXc- 



(37) 



r^c^c^r 



Using these relations we find for the optimal permuta- 
tion that dL(p2) - °-^^'^''' A'^A'_„,, \v 12 
-S^dtma-Xmn \X., 
we need to prove that 



-87di[2(l - A')A'max™„ \Xmni j - 
^L[p]. To prove \dL{p)\ < \dL{p2)\ 



At- Ap LAr 



/ ^ y\r^^c I 
r.c 



< 



/ ^ ^r^c 



max \X„ 

mn 



(38) 



This is trivially true since all the A^'s are positive. All 
that remains is to bound the maxm„ |^mnP in any 
unbiased-basis. 

A general unitary that transforms the basis \k) to an 



(D-ir 



E\X P I (33) unbiased-basis is T\n) = ^''^j ^^cxp{i<j)i!^')\k). It is 



r^p,c^p 



possible to rewrite iX^^nP as |X„ 



|(m|rtj^T|n)| = 



so that 



niax|X„j„p < max — -j 



J2 e*(^i"'-<"")fc 



< 



£12 






(39) 



For even D Eq. dMl) evaluates to D'^/16; for odd D it 
evaluates to D^/lG - 1/8 + 1/(16L'2). For large D (say 
D > 5), D^/16 is a good approximation for both even 
and odd D. To find the lower bound on the decrease in 
impurity it is important to remember that two matrix el- 
ements contribute to the sum: niaxm.„ \Xmn\'^^m'^n, and 
maxm^n |^nmpA„Am. For large D the impurity under the 
two eigenvalue distribution is then 



dL:> < 



-8jdt—2XoXi 
16 



-^dt^L2{t). (40) 



The dependance for this matrix element is only on the 
dimension of the system. Thus the speed-up upper bound 
for any unbiased-basis feedback is 



S2< 



2 ' 



(41) 



for D > 1 and i > 7" 



V. A REGISTER OF QUBITS 

We now generalize UBB feedback protocols to the case 
of a register of n qubits, where each qubit is indepen- 
dently and continuously measured. Instead of one ob- 
servable X, we now have n, given by X^''^ = J*^^) (g) /f^^) (^ 
. . . cri''-' . . . (g) /("' , where r labels the rth qubit. The SME 
describing such a measurment is 

dp = ^2Kdtv\x^'^A p + V2^dW^''^n [x('')j p. (42) 



The combined state of the n qubits exists in a I? = 2"^' 
dimensional Hilbert space. 



A. Purification from Measurement Alone 

In this section we will not analyse the no-feeback case 
to the same level of rigor as we did in section III Al but 
rather, we rely upon the intuition gained from that anal- 
ysis. For simplicity, consider first a two-qubit register 
with uncorrelated qubits. The state of the system is 
p ~ Pi® P2- The impurity for such a state is 

i(2) = 1 _ Tr [(pi ® p2f] = 1 - Tr [p?] Tr [pi] , (43) 

where the superscript "(2)" on L signifies the number of 
qubits in the register. For very pure states it is natural 



to parameterize the eigenvalues of the rth qubit, pr, as 
{Ag = 1 — Ar, A^^ = Ar} where the convention that A^ <C 
1 still holds. To first order in A we have 

L(2) = 1 - (1 - 2Ai)(l - 2A2) - 2(A2 + Al). (44) 

Because the qubits are initially uncorrelated, and are in- 
dependently monitored, it is reasonable to assume that 
Al and A2 are of the same order, so L'^^) ^ 2(2A). When 
A < 1 then L = 2(1 - A)A « 2A so L^^) « 2L. At long 
times, for n qubits we expect that 



{L 



{n)^ 



LT 



l{L 



(1)^ 



LT- 



(45) 



We can verify this using the analytic expression for qubit 
impurity, Eq. ([7|) . Recalling that the purity may be writ- 
ten as P = 1 - L, we have (p(")) = (p(i))", so 



(L(")(t)) 



-iKt 



-RVm^ji 



87rt J-00 cosh(v8Ki?) 



(46) 



Using the long-time expression for a qubit, Eq. ([5]), we 
obtain a long-time analytical expression for L in the ab- 
sence of feedback for a register of n qubits 



(L(")(i)) 



LT 



1 



1 



nne 



-4Kt 






-4Kt 



V647r7t ^ 



(47) 



To compare this expression to the qubit expression, 
Eq. ([U, n is set to 1 and k ~ j-f. 



B. Lower bound for purification using UBB 
feedbaclc on a register of qubits 

The change in impurity for a register of qubits in an 
unbiased basis is 



n 

dL = -SKffi^Tr [xWp^^'' 



(48) 



n (-D-1) 

= -SKdtY, Y. \^S\'^^^J- (49) 

r— 1 iJ—0 

Using pf to calculate the lower bound for a UBB pro- 
tocol in a register appears to be difficult. Instead we 
will derive the lower bound using the same method pre- 
sented in Appendix |B] To do this we will again in- 
troduce the feedback in the Heisenberg picture so that 
Xir.Tn) ^ p^TX^^'^T^Pl. As before the T's are condi- 
tional unitaries that introduce the unbiasedness (between 
p and X'^'')) and the Pm's also retain their meaning as 
permutations. For a register being measured in a basis 
that is unbiased with respect to the logical basis, with a 
randomly changing permutation, the change in impurity 
is 



dL 



n D! 

EE 



-SKdtTv 



x(''.")px('''"V 



where i) = 2" — 1 . After performing a similar procedure 
to the one found in Appendix [B] we find 



dL = -SKdtnD{D-2)\L{t). 



(50) 



Thus we find the impurity of the state undergoing feed- 
back of the above form decreases as 



i(")(t) 



-8Knt/{D- 



i)L(O). 



The asymptotic speed-up factor is 



D- 1 



(51) 



(52) 



For n = 1 (i.e. a qubit) the speed-up factor 5' is S* = 2, 
which agrees with the resuh from Ref. [30|. When n = 2, 
the speed-up is S* = 4/3. This is comparable to the speed- 
up found for the locally optimal rapid measurement (RM) 
protocol in |4l| where the predicted speed-up in the long- 
time limit was S-rm ~ 1-4. (The feedback in the rapid 
measurement protocol permutes the eigenvalues of p in 
the logical basis to decrease a different measure of impu- 
rity. Thus at all times the state eigenbasis and the mea- 
surement basis commute.) Unfortunately when n = 3, 
S = 6/7, which is a slow-down. This slow-down trend 
continues for all ri > 3, and for large n the slow-down is 
~ n2-"+i. 

One can interpret Eq. ()52p as a lower bound on an 
UBB algorithm for a register of qubits in much the same 
way as we did in Section fill Bl Equation ([5^ may also 
be interpreted as an all-permutation deterministic purifi- 
cation protocol, although it is only useful in a two qubit 
register. 



C. Upper bound for purification using UBB 
feedback on a register of qubits 



As before it is possible to rewrite |Xmn| as |Xmri| 



|(m|Ttx('')r|ri)| so that 



maxlX^J < max — 



ran 



D-1 



y^ e*'^''(— 1)-^^'^''"^ 



k=0 



, (53) 



where /(/c, r) is a function that appropriately determines 
the sign of the diagonal elements of AT^'"^. Thus 



(54) 



and 



dL^2 < -8Kdt'^L{t) = -%KdtnL{t). 



At long times the speed-up upper bound for any UBB 
feedback in a register is 



By substituting in n = 1 wc regain the result of Ref. [30| . 
At present it is unclear if the bound in Eq. ([55]) is tight. 
One reason why we expect this bound not to be tight 
comes from the upper bound on the rapid measurement 
protocol of Ref. [41[ . The upper bound was shown to be 
n and numerical results indicated it was S-rm = 0.718n, 
which is significantly less than the upper bound. 



VI. DISCUSSION 

For a qudit we have bounded the speed-up in purifi- 
cation for any unbiased basis feedback by |(D -I- 1) < 
S < D^ /2. Prior to this work only the lower bound was 
known for the qudit. In Ref. [l9| it was shown that by 
applying a unitary, chosen at random from the unitary 
group il(D), to a qudit during the measurement process 
(random feedback) one could attain 5* = |D asymptot- 
ically. This suggests UBB feedback is unnecessary for 
a speed-up 0{D), although it remains an open question 
whether the actual speed-up for UBB protocols is 0{D) 
or 0{D^). 

We also examined UBB feedback applied to a register 
of n qubits. We have shown that the speed-up in purifi- 
cation is bounded by 2n/(2" — 1) < 5 < 2n. Prior to this 
work there were no results for purification of a register 
of qubits using UBB feedback. 

Unlike the rapid measurement (RM) protocols of [41| 
the UBB protocols presented in this paper are not lo- 
cally optimized in time for decreasing dL. Additionally 
it is not clear if UBB protocols allow one to obtain in- 
formation about the preparation procedure, unlike RM 
protocols which may be used for state estimation. 

In the future we plan to investigate locally optimal 
protocols for the impurity (and other measures of mixed- 
ness), and examine rapid purification using a particular 
unbiased basis. Also, very recently Ruskov et at. [4g 
have considered monitoring a single qubit in three un- 
biased bases simultaneously. The relation between their 
work and ours remains to be explored. 
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Appendix A: Proof that Xu — 

For a traceless operator AT, which is diagonal in the 
computational basis, we can find a transformation T that 
will make it unbiased with respect to the computational 
basis. The transformation has the form X = TXT^ 



S2 < 2n. 



(55) 



Xii — 



D 



'.T, 



Now we wish to examine the diagonal elements. Being 
mindful of the fact that Tj,, = T* we have 



rs sr 



D D D 

Xu = 2^ TijXjkTi^. = 2^ T^jXjkT*^, = 2^ T.,jT*,^Xjk- 

j,k=l j,k=l j,k=l 

Using the fact that X is diagonal gives the expression 



xu = Y.T^,^x,,. (Ai) 

Unbiasedness means that TijT*- = 1/D so that 



1 ^ 






(A2) 



Recalling that X is traceless, we thus have Xa = 0. 



E,(a|(l.7)01)l&))sothat 

^ - E E A,„(.)(X*|j)0-|P^pP„|X*) (B4) 

rn— 1 i j' — 1 

D\ D 
= E E -^>n(OA„(,)(Xi|j>0-|Xi) (B5) 

rn— 1 iJ — 1 

D\ D 

= E E K^i^)K.U)\{X^\J)?■ (B6) 

rn— 1 i,^ — 1 

Next we wish to sum over the m label, which represents 
permutations of the basis. First we define 



D! 



Q, = |(Xz|j)|2^A„(,)A„ 



0)' 



so that S = ^- Cij. Suppose i — j, then 



D! 



cu = \{m)\'J2>^l 



,(i)- 



(B7) 



(B8) 



Appendix B: Proof of Eq. ([26ll 

Denote eigenvalues and eigenvectors of p in the basis 
which diagonalizes p as p|i) = Xi\i). Define X as a Her- 
mitian operator that is diagonal in this eigenbasis, and 
X = TXT^ as a transformed version of X that is unbi- 
ased with respect to this eigenbasis. Now define P„j, for 
m ~ 1, . . . ,Dl, to be the D\ operators that give each of 
the possible permutations of the eigenbasis \i): for ex- 
ample, Pm\i) = l'Ti(i)) labels a particular permutation of 
the basis \i), while P^K) — |?ti~^(*)) is an inverse per- 
mutation to m. Note that P^Pm|z) = PmP^|z) = \i). 
Starting with Eqn. [26] one may write 



D\ 

S = E Tr [P™XPipP,„XP,t„p] 

m— 1 
D\ D 

= E Y.^^xp^^pp„,xp},,pp„. 



(Bl) 



m— 1 z— 1 



Consider the action of the permutation operators and 
the state on the basis \i) 

P}nPP^n\i) = Pl,p\m{i)) = A,„(,)P,|jm(i)) = \m(€)\i) ■ 

(B2) 
Using the notation X\i) = \Xi), we thus have 

D\ D 

S ^ Y.T.^"^i^}^^'\P^nPPn.m. (B3) 

m— 1 z— 1 

Now we use the completeness of the basis {{a\b) = 



There are {D — 1)! permutations that take i -^ p for a 
fixed pe [1,2..., D], so 

D 

C,. = |(X^|^)p(i?-l)!EA^-|(^*N)P(^-l)!Tr[p^]. 
p=i 

(B9) 
Now consider the case where i y^ j. Take any pair 
q ^ p and look for permutations that take i ^ p, j ^ q. 
There are {D — 2)1 of these permutations: 

D 

a, = \{Xi\j)\\D-2y.Y,\K (BIO) 

Now examine J2pj,q ^p'^g term: 

(E^p)' = (Eap)(Ea,) 

P P Q 

D D 

p pi^q 

D 

(Tr[p]f = Tr[p2]+^AA. (Bll) 

p¥=q 

As p is normalized we find X]p=£o W = 1 — Tr [p^] . 
The simplified expression for Eq. (|B1[) . so far, is 



s = E^^ + E^'^ 

= Y^{D-\)nr{p\\{X^\i)\^ 

i 

+ E(^-2)!(l-Tr[p2])|(X^b•)p. (B12) 



Now we simplify the expression J2i^i l(^*b)P- Us- 
ing the Parseval relation for a vector ^, (^|^) = 



-D /v„-i„-\2_T. r^i2 



^D 



E,"il(*|j>P,wefind 

D 



D 



3=1 



D 



J=» 



{X^\X^) = Y.U^\m^ + {X^\i)^ (B13) 



^D 



i-e- Ej^^i l(^i|i>l = (^i|^») - (Xi|i)2. The total ex- 
pression is now 

D 

S = ^(Z?-l)!Tr[p]|(X^|^)p 

i 

+ {D 2)!(1 - Tr [p^])f2{{X^\X^) {X^\^)^). 

(B14) 
Massaging the Ei=i(^*l^*) term gives 

D D 

Y^mXi) = ^(z|XtX|z)=Tr[x2]. 

RecaU that X = TXT^, thus Tr [X"^ X] = 
Tr [(TXrt)trXTt] = Tr [X^] . 

Now we examine the Ei=i("^*K)^ term. The eigen- 
vectors and eigenvalues of X are Xj^I/^) = Xb\'^b) and 
{i\^b) = ^l^fD by way of their unbiasedness. By insert- 
ing the identity we have 



{Xi\i) = Y,{Xi\^b){^ 



bW 



b 

^Y.^b(i\^b) 

b 

^;^a.b=Tr[X]/Z?, (B15) 



so Et.i(^*N)^ = Tr [X] /D^. The total expression is 



wn2 



5 = (Z?-l)!Tr[p]Tr[X] /i?2 

+ p - 2)!(1 - Tr [p']){Ti [^2] - Tr [X]" /D''). 

(B16) 

Finally, recalling from Appendix [^ that X is traceless, 
we have 

S ^ {D^ 2)!Tr [X^] (1 - Tr [p^] ) (B17) 



as required. 



Appendix C: Proof of Eq. dM] 



D \v 12 _ Y--D \v 12 



Recall that Er^o l^-oP = Er=o l^-oP as |A,,| = 0. 
So, 

(fc|Al"A|fc) = (fc|(t/J,[/^)1"[/J^t/^|fc) == (fc|[/jft/^|fc). 

Now insert the identity 

{k\X^X\k) = Y.{k\U\l){l\Jl\m){m\U^\k) 

1.71% 

= Y.{k\U\l){l\Jl\m){m\U^\k). (CI) 

l.ra 

The matrix J^ is diagonal, so 

{k\x^x\k) = 5](/|j2|0(fc|[/|0(/|[/^|fc) 

= Y.^i\Jl\i)mu\i)? 



D 



12 
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